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ABSTRACT 

Polarization  diversity  in  radar  is  useful  to  detect  targets, 
particularly  when  Doppler  discrimination  is  not  possible. 
We  introduce  a  new  polarimetric  radar  model  that  includes 
the  realistic  dependence  of  the  clutter  reflections  on  the 
transmitted  signal.  We  then  develop  a  statistical  polarimet¬ 
ric  detection  test,  robust  to  heavy  inhomogeneous  clutter. 
We  prove  that  the  detector’s  false-alarm  rate  is  invariant 
to  the  space  and  time  variabilities  of  the  clutter;  hence,  it 
has  the  constant  false-alarm  (CFAR)  property,  while  still 
maintaining  a  good  probability  of  detection.  We  demon¬ 
strate  the  improved  performance  of  the  proposed  detector 
in  comparison  with  existing  detectors  using  numerical  ex¬ 
amples  and  real  data. 

1.  INTRODUCTION 

The  detection  of  static  or  slowly  moving  targets  in  heavy 
cluttered  environments  is  considered  a  challenging  prob¬ 
lem,  mainly  because  it  is  not  possible  to  discriminate  the 
target  from  the  clutter  using  the  Doppler  effect.  Polariza¬ 
tion  diversity  provides  additional  information  that  enhances 
the  detection  of  targets,  particularly  under  the  conditions  de¬ 
scribed  above. 

Earlier  work  in  the  field  of  polarimetric  detection  has 
addressed  the  design  of  detectors  under  the  assumption  that 
training  data  is  available  [l]-[3].  The  performance  of  these 
detectors  can  be  severely  degraded  in  the  presence  of  in¬ 
homogeneous  and  non-stationary  clutter.  To  overcome  this 
problem,  the  application  of  compound-Gaussian  distribu¬ 
tions  for  modeling  the  global  behavior  of  nonhomogeneous 
clutter  has  been  proposed.  However,  the  use  of  non-Gaussian 
models  increases  the  difficulty  of  developing  efficient  detec¬ 
tion  algorithms.  For  instance,  the  detection  statistic  derived 
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from  the  scheme  presented  in  [4]  has  a  closed-form  expres¬ 
sion  only  for  the  special  case  of  two  polarimetric  channels, 
and  the  detector  proposed  in  [5]  does  not  support  the  con¬ 
stant  false-alarm  rate  (CFAR)  property  with  respect  to  the 
clutter  covariance. 

In  this  paper,  we  develop  a  polarimetric  detector  whose 
decision-making  capability  is  based  on  the  data  collected 
only  from  the  range  cell  under  test,  without  resorting  to  sec¬ 
ondary  data  or  prior  knowledge  of  the  target  and  clutter.  In 
order  to  derive  the  detector,  we  first  introduce  a  new  polari¬ 
metric  radar  model  that  states  the  dependence  of  the  clut¬ 
ter  reflections  on  the  transmitted  signal.  We  prove  that  it 
is  robust  against  inhomogeneous  and  non-stationary  clutter; 
i.e.  the  detector’s  false-alarm  rate  is  insensitive  to  space  and 
time  variations  of  the  clutter,  while  still  maintaining  a  good 
probability  of  detection.  Moreover,  this  test  is  CFAR.  In 
that  sense,  this  detector  resembles  that  of  [6];  however,  our 
work  is  distinctive  in  several  ways:  (i)  The  detector  in  [6] 
is  developed  under  the  assumption  that  the  target  is  present 
only  in  one  of  the  multiple  scans  of  the  radar  dwell;  thus, 
it  will  not  detect  a  static  target.  Our  method  is  specifically 
aimed  at  detecting  static  or  slow  targets,  (ii)  Our  detector 
has  been  extended  to  consider  any  arbitrary  transmit  polar¬ 
ization.  (Hi)  Our  detector  can  also  be  applied  to  advanced 
polarimetric  receivers,  such  as  tripoles  and  vector-sensor 
antennas,  in  addition  to  conventional  vertical  and  horizontal 
(or  left  and  right  circular)  polarized  arrays  of  antennas,  (iv) 
Our  detector  is  tested  with  real  radar  data.  An  additional 
advantage  of  our  detection  algorithm  is  that  it  is  computa¬ 
tionally  less  intense  because  it  does  not  require  processing 
of  training  data. 

2.  POLARIMETRIC  RADAR  MODEL 

We  consider  a  mono-static  radar  capable  of  transmitting  wave¬ 
forms  with  any  arbitrary  polarization  on  a  pulse-by-pulse 
basis.  In  addition,  we  assume  that  the  radar  system  illumi¬ 
nates  a  point  target.  The  recorded  data  from  the  range  cell 
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under  test  consist  not  only  of  the  target  echoes  but  also  of 
the  undesired  reflections  from  the  target  environment  (clut¬ 
ter). 

2.1.  Polarimetric  Data 

The  output  of  a  diversely  polarized  array  of  Q  sensors  re¬ 
ceiving  the  echoes  from  the  cell  under  test  can  be  expressed 
as 


•  The  vector  e(t)  represents  the  thermal  noise. 

Equation  (1)  can  be  written  as 

y(t)=s(t)B£(/j,  +  x)  +  e(t),  (5) 

where  the  scattering  coefficient  vectors  of  the  target  and 
clutter  are  y  =  [4h>  4v>  shvF  and  x  =  [s£h,  s=v,  s£v]T, 
respectively,  which  have  dimension  P  =  3.  The  polariza¬ 
tion  matrix  £  is 


y(t)  =  B  (S'*  +  Sc)  m  +  e(i),  t=l,...,N,  (1) 

where 

•  The  Q  x  1  vector  y(t)  is  the  complex  envelope  of  the 
measurements. 


•  The  Q  x  2  matrix  B  is  the  response  of  the  diversely 
polarized  sensor  array.  If  the  receiver  array  is  a  vector 
sensor  [7],  the  array  response  is  given  by 

—  sin  <f> 
cos  <p 

B=  ° 

—  cos  <f>  sin  ip 

—  sin  <p  sin  ip 

cos  ip 

where  <p  and  ip  are  the  azimuth  and  elevation  angles 
of  the  cell  under  test,  respectively.  For  a  conventional 
polarized  radar  measuring  the  horizontal  and  vertical 
components  of  the  electric  field  and  assuming  these 
two  sensors  are  orthogonal  to  the  direction  that  points 
toward  the  cell  under  test,  the  array  response  matrix 
is  B  =  I2,  where  I2  is  the  2x2  identity  matrix. 

•  The  complex  scattering  matrix  S  represents  the  po¬ 
larization  change  of  the  transmitted  signal  upon  its 
reflection  on  the  target  or  clutter: 


cos  <p  sin  ip 
-  sin  <p  sin  ip 
cos  ip 
sin  <p 
—  cos  <p 
0 


(2) 


S  = 


Sll 

S21 


(3) 


where  the  variables  Sn  and  S22  are  co-polar  scatter¬ 
ing  coefficients  and  S12  and  S21  are  cross-polar  coef¬ 
ficients.  For  the  mono-static  radar  case,  S12  =  S21. 
The  superscripts  *  and  c  refer  to  the  target  and  clutter. 

•  The  vector  £(f)  is  the  narrowband  transmitted  polar¬ 
ized  signal  which  can  be  represented  by 


m 


cos  a  sin  a 

cos  (3 

—  sin  a  cos  a 

j  sin  / 3 

(4) 

where  £1  and  £2  are  the  signal  components  on  the  po¬ 
larization  basis  of  the  transmitter,  a  and  (3  are  the  ori¬ 
entation  and  ellipticity  angles,  respectively,  and  s(t) 
is  the  complex  envelope  of  the  transmitted  signal. 


£  = 


&  0  6 

0  £2  Ci 


(6) 


The  time  samples  can  be  stacked  in  one  vector  of  dimension 
NQ  x  1: 

y  =  (s  <g>  B%)  {y  +  x)  +  e,  (7) 

where  s  =  [s(l), . . . ,  s(N)]T  and  <g>  is  the  Kronecker  prod¬ 
uct.  Piling  together  the  data  corresponding  to  pulses  of  dif¬ 
ferent  polarization  yields 


y  =  Ay  +  Ax  +  e, 


(8) 


where  A  is  a  M  x  P  (M  =  KNQ)  complex  matrix  that 
represents  the  system  response: 


s  <8>  Bt;  1 
s  (g)  B^k 


(9) 


and  Cfc  is  the  polarization  matrix  of  each  diversely  polarized 
pulse  (k  =  1, . . . ,  K). 


2.2.  Statistical  Model 

We  assume  that  the  target  is  a  small  man-made  object;  hence, 
y  is  a  deterministic  vector.  On  the  other  hand,  the  clutter  in 
the  range  cell  under  test  can  be  considered  as  a  large  collec¬ 
tion  of  point  scatterers  producing  incoherent  reflections  of 
the  radar  signal.  Then,  x  is  a  zero-mean  complex  Gaussian 
random  vector  with  covariance  matrix  E.  The  noise  e  is 
a  zero-mean  complex  Gaussian  random  vector  with  covari¬ 
ance  matrix  gIm  ,  where  Im  is  the  M  x  M  identity  matrix. 
In  addition,  we  assume  that  the  clutter  reflections  and  the 
thermal  noise  are  statistically  independent. 

The  radar  dwell  consists  of  a  series  of  pulses  that  can 
be  seen  as  “snapshots”  of  the  range  cell  under  test.  If  the 
pulse  duration  is  short  with  respect  to  the  dynamic  of  the 
target  and  its  environment,  it  is  reasonable  to  assume  that 
their  scattering  coefficients  are  constant  during  each  pulse. 
However,  from  pulse  to  pulse,  we  consider  the  clutter  scat¬ 
tering  coefficients  as  independent  realizations  of  the  same 
random  process.  Then,  the  distribution  of  each  snapshot  is 

yd  ~  CN(Ay,  AYiAh  +  <jIm),  d=l,...,D,  (10) 

where  D  is  the  total  number  of  snapshots  or  pulses  in  the 
radar  dwell  and  H  denotes  the  conjugate  transpose  operator. 


2.3.  Known  and  Unknown  Model  Parameters 


The  MLE  of  E  is 


We  assume  the  system  response  matrix  A  is  known,  since 
we  consider  that  the  receiver  antenna  array  has  been  cali¬ 
brated.  In  addition,  we  assume  that  the  power  of  the  ther¬ 
mal  noise  er  is  known,  because  it  can  be  easily  estimated 
from  the  recorded  data  when  no  signal  has  been  transmit¬ 
ted.  However,  we  suppose  that  we  have  no  prior  knowledge 
about  the  target  and  the  clutter,  nor  do  we  count  on  a  sec¬ 
ondary  data  set  for  estimating  the  statistical  properties  of 
the  clutter.  Hence,  the  vector  //,  and  the  matrix  E  are  the 
unknown  parameters  of  the  statistical  data  model  (10). 


3.  DETECTION  TEST 


We  aim  to  decide  whether  a  target  is  present  or  not  in  the 
range  cell  under  test,  based  on  the  recorded  data.  Then,  the 
decision  problem  consists  of  choosing  between  two  possible 
hypotheses:  the  null  hypothesis  Ho  (target-free  hypothesis) 
or  the  alternative  hypothesis  Hi  (target-present  hypothesis) 


f  W0:  H  =  0,  E 

\Hi:  y  ±  0,  E 


(11) 


where  the  matrix  E  is  considered  as  a  nuisance  parameter. 
Next,  we  derive  the  generalized  likelihood  ratio  (GLR)  test 
[8]  and  study  its  performance.  We  omit  derivation  details 
in  this  paper  due  to  space  limitation;  refer  to  [9]  for  an  ex¬ 
tended  version  of  this  paper. 


3.1.  GLR  Test 

The  GLR  test  decides  Hi  if 


£glr 


/o(3/i,---,3/d;£o) 


(12) 


where  /o  and  fi  are  the  likelihood  functions  under  Ho  and 
Hi,  So  and  Ei  arc  the  maximum  likelihood  estimates  (MLE) 
of  E  under  Ho  and  Hi,  /h  is  the  MLE  of  //  under  Hi,  and  7 
is  the  detection  threshold.  For  simplicity  of  notation  we  will 
omit  references  to  the  data  in  the  arguments  of  the  functions 
fo  and  fi  in  the  rest  of  the  paper. 

We  start  the  derivation  of  the  GLR  test  by  determining 
the  likelihood  functions  and  MLEs  of  the  unknown  param¬ 
eters.  Under  hypothesis  Ho ,  it  is  assumed  that  fi  -  0;  then 


ln/0(E)  =  -D  [Mhnr  +  ln\C\  -(- tr(C_1S0)]  ,  (13) 


E0  =  A+S0A+H  -  a{AHA)-i,  (15) 

where  A+  =  (AH A)~x AH  is  the  generalized  matrix  in¬ 
verse.  Under  hypothesis  H  \ ,  the  likelihood  function  is 


ln/i(/x,E)  =  —D  |^Mln7r  +  In  \C\  +  tr(C,_1C'1) 
where 


Ci  =  4  (yd  ~  A p)  (yd  ~  ap)h  • 


D 


D 


d—1 


The  MLE  of  the  unknown  parameters  are 

Mi  =  A+y, 

Ei  =  A+SiA+H  -a(AHA)~1, 
where  y  is  the  sample  mean  vector 


D 


y=pJ2ydi 


d=  1 


and  Si  is  the  sample  covariance  matrix 


D 


si  =  -p  (yd  -  y)  (yd  -  v)H  ■ 


(16) 


(17) 


(18) 

(19) 


(20) 


(21) 


d=  1 


Then,  substituting  the  MLEs  and  likelihood  functions  in 
(12),  and  after  some  mathematical  simplifications,  the  GLR 
test  can  be  written  as 


Cglr 


1  +  yHA  (Ah SiA)  1AHy 


(22) 


Since  (22)  is  a  monotonically  increasing  function  of  the  sec¬ 
ond  term  inside  the  brackets,  an  equivalent  detection  test 
statistic  can  be  defined  as 

Tglr  =  yH A  (AhSiA)~1  AHy.  (23) 


3.2.  Detection  Performance 


Applying  Corollary  5.2.1  from  [10],  it  is  straightforward  to 
verify  that  the  detection  statistic  (23)  is  distributed  as  fol¬ 
lows 


Tglr 


D-P 

P 


p2P,2(D-P)  under  Ho 

•77,2p,2(d-p)(A)  under  Hi 


where  C  =  AT,AH  +  ctIm  is  the  theoretical  covariance 
matrix  of  the  data,  defined  in  (10),  and  ,S'o  is  the  sample 
covariance  matrix 


d= 1 


(14) 


where  Pvx,v2  denotes  an  F  distribution  with  ui  and  V2  de¬ 
grees  of  freedom,  and  F' Vl.V2{X)  denotes  a  non-central  F 
distribution  with  Vi  and  V2  degrees  of  freedom  and  non-centrality 
parameter  A.  The  non-centrality  parameter  is  given  by 


A  =  2  DfiH 


H  +  a(AHA)  1 


M 


(25) 


Figure  1:  Detection  performance  of  the  GLR  test  statistic 
Tglr  as  a  function  of  the  target-to-clutter  ratio,  for  different 
values  of  Pfa- 

Thus,  the  detection  performance  becomes 

ft  =  Qr>  2P,2(D-P)(A)(7)»  (26) 

where  1\)  is  the  probability  of  detection,  PpA  is  the  proba¬ 
bility  of  false  alarm,  Q  is  the  right-tail  probability  function 
[8,  Chap.  2]  and  7  is  the  detection  threshold  for  the  required 
probability  of  false  alarm.  In  particular,  note  that  the  ex¬ 
pression  for  Pfa  does  not  depend  on  the  covariance  of  clut¬ 
ter  and  thermal  noise,  nor  on  the  transmitted  signal;  hence, 
equation  (23)  is  a  CFAR  test. 

We  validated  the  analytical  performance  of  the  devel¬ 
oped  test  by  computer  simulations.  We  performed  Monte 
Carlo  simulations  based  on  100/ Pfa  independent  trails.  Fig¬ 
ure  1  shows  curves  of  Pd  as  a  function  of  target-to-clutter 
ratio  (TCR)  for  different  values  of  Pfa  and  a  clutter-to- 
noise  ratio  (CNR)  of  lOdB.  The  definitions  of  TCR  and 
CNR  are  given  in  the  Appendix.  The  results  depicted  in  the 
figure  consist  of  the  average  of  100  different  cases.  For  each 
case,  the  entries  of  //,,  E,  and  A  are  independent  realization 
of  a  random  variable  with  distribution  CJ\f{ 0, 1);  then  those 
vectors  and  matrices  are  scaled  to  fulfill  the  required  TCR 
and  CNR. 

3.3.  Detector  Performance  in  Inhomogeneous  Clutter 

A  significant  characteristic  of  the  detector  presented  in  this 
paper  is  that  the  test  statistic  (23)  is  computed  using  only 
data  from  the  cell  under  test.  No  assumption  is  required 
about  the  range  homogeneity  of  the  clutter.  Thus,  the  detec¬ 
tor  is  robust  against  spatially  fluctuating  clutter. 

To  support  this  assertion,  the  Pfa  is  computed  using 
5  •  105  runs  of  Monte  Carlo  simulations  in  the  presence  of 


inhomogeneous  clutter.  The  clutter  echoes  are  assumed  to 
follow  a  compound-Gaussian  distribution  with  covariance 
matrix  rE.  The  texture  r  has  a  generalized  Gamma  proba¬ 
bility  density  function  [4] : 

p(T)  =  f^G)l>"1“p(-r)’  <27) 

where  v  is  the  order  parameter  and  S  is  the  average  power. 
Note  that  v  =  00  corresponds  to  clutter  echoes  with  Gaus¬ 
sian  distribution,  but  when  v  decreases,  they  deviate  from 
Gaussian.  Additional  information  about  the  simulation  setup 
is  as  follows.  The  transmitted  signals  are  rectangular  pulses, 
which  are  transmitted  in  a  sequence  of  alternating  vertical 
(V)  and  horizontal  (H)  polarization  (K  =  2).  The  radar 
dwell  consists  of  10  repetitions  of  this  polarimetric  sequence 
( D  =  10).  On  the  receiver  side,  samples  of  the  V  and  H 
electric  field  are  recorded  ( Q  =  2)  at  a  rate  of  one  sam¬ 
ple  per  pulse  (N  =  1).  Since  the  simulations  are  intended 
to  compute  the  Pfa.  no  target  is  considered.  The  average 
power  of  the  texture  is  S  =  50  and  the  speckle  component 
of  the  covariance  is  E  =  I3. 

For  comparison  purposes,  we  also  compute  the  perfor¬ 
mance  of  other  detectors  that  use  secondary  data: 

•  Polarization-space-time  GLR  (PST-GLR):  This  detec¬ 
tor  was  derived  under  the  assumption  of  homogeneous 
Gaussian  clutter  [1], 

•  Texture  free  GLR  (TF-GLR):  This  detector  was  for¬ 
mulated  assuming  that  the  clutter  follows  a  compound 
Gaussian  distribution  [4], 

•  Hybrid  U  test;  This  detector  was  developed  for  a  sin¬ 
gle  polarimetric  channel  assuming  that  the  clutter  fol¬ 
lows  a  Weibull  distribution  [11], 

In  Fig.  2,  the  Pfa  of  our  detection  test  Tglr.  as  well  as 
of  the  PST-GLR,  TF-GLR,  and  Hybrid  U  tests,  is  plotted  as 
a  function  of  the  order  parameter  of  the  clutter  texture  dis¬ 
tribution.  For  the  latter  three  tests,  we  consider  80  adjacent 
range  cells  for  generating  the  secondary  data.  We  observe 
that  the  performance  of  our  detector  Tglr  and  the  TF-GLR 
test  remains  constant  at  the  designed  false-alarm  rate.  How¬ 
ever,  there  is  a  discrepancy  between  the  design  and  the  ac¬ 
tual  clutter  that  increases  the  probability  of  false  alarm  for 
the  PST-GLR  test  as  the  clutter  process  departs  from  being 
Gaussian.  The  Hybrid  U  test  behaves  similar  to  the  PST- 
GLR;  but  it  has  larger  Pfa  due  to  the  mismatch  between  the 
model  used  to  generate  the  clutter  data  and  the  one  assumed 
by  the  algorithm. 

3.4.  Detector  Performance  Using  Real  Data 

To  evaluate  the  performance  of  the  polarimetric  detectors 
in  a  real  application  situation,  we  use  data  collected  at  the 


Figure  2:  Probability  of  false  alarm  as  a  function  of  the  or¬ 
der  parameter  v  of  the  clutter  texture  distribution. 

Osborne  Head  Gunnery  Range  (OHGR),  Dartmouth,  Nova 
Scotia,  Canada,  with  the  McMaster  University  IPIX  radar 
[12],  Specifically,  we  use  the  data  recorded  on  November 
11,  1993.  Fig.  3  shows  the  magnitude  of  the  echoes  from 
a  small  target  in  inhomogeneous  sea  clutter  for  the  four  po- 
larimetric  channels:  HH,  VH,  VV,  and  HV. 

Fig.  4  shows  detection  maps  for  the  Tqlr,  PST-GLR, 
TF-GLR,  and  Hybrid  U  tests.  The  detection  threshold  is 
again  set  for  PpA  =  10~3,  and  other  parameters  are  also 
as  in  the  former  subsection.  For  the  latter  three  tests,  a 
guard  region  of  two  range  cells  on  either  side  of  the  cell 
under  test  is  used  to  avoid  the  leakage  of  the  target  energy 
into  the  secondary  data.  We  observe  that  our  detection  test 
finds  the  true  position  of  the  target  with  only  a  few  false 
detections  despite  the  strong  clutter  echoes,  which  can  be 
appreciated  from  Fig.  3.  On  the  other  hand,  the  PST-GLR 
has  high  false-alarm  rate  due  to  the  presence  of  inhomoge¬ 
neous  clutter,  as  was  predicted  in  the  previous  subsection. 
The  Hybrid  U  and  TF-GLR  tests  reduce  the  false-alarm  rate 
with  respect  to  the  PST-GLR  test;  however,  it  still  remains 
considerably  high. 

One  possible  way  of  summarizing  the  performance  of 
these  detection  tests  in  relation  to  real  data  is  to  plot  the  em¬ 
pirical  receiver  operating  characteristic  (ROC)  curves.  The 
empirical  ROC  is  obtained  by  dividing  the  detection  map 
in  two  regions.  One  region  corresponds  to  the  range  where 
the  target  is  located  and  the  other  region  corresponds  to  the 
range  where  there  is  only  clutter.  For  different  values  of 
the  detection  threshold  7,  the  empirical  Pd  and  Pfa  are 
computed  by  counting  the  pixels  on  each  respective  region. 
Fig.  5  shows  that  the  empirical  ROC  curve  of  the  Pglr  test 
is  on  the  top,  outperforming  the  other  algorithms. 

We  acknowledge  that  comparing  the  Hybrid  U  and  TF- 
GLR  tests  with  the  other  two  is  not  strictly  fair,  since  these 
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Figure  3:  Magnitude  in  range  and  time  domain  of  the  IPIX 
radar  dataset  stare  1,  collected  on  November  11,  1993.  The 
target  location  is  indicated  by  the  markers  “>”  and  “<”. 
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Figure  4:  Detection  map  in  range  and  time  domain  for  the 
IPIX  radar  dataset  starel,  collected  on  November  11,  1993. 
Black  pixels  indicate  that  the  detection  statistic  is  larger  than 
the  threshold.  The  target  location  is  indicated  by  the  mark¬ 
ers  “>”  and 


Probability  of  False  Alarm 

Figure  5:  Empirical  receiver  operating  characteristic  (ROC) 
curves. 


tests  were  developed  for  one  and  two  polarimetric  channels, 
respectively.  (We  have  used  HH  and  VV  for  the  previous 
results  of  the  TF-GLR  test.)  Our  goal  is  not  to  extend  these 
algorithms  but  to  show  that  detection  tests  using  secondary 
data  may  have  larger  false  false-alarm  than  their  designed 
value,  in  the  case  of  inhomogeneous  clutter. 

4.  CONCLUSIONS 

In  this  paper,  we  addressed  the  problem  of  detecting  small 
and  slowly  moving  targets  in  heavy  inhomogeneous  clutter 
by  exploiting  polarization  diversity.  We  first  introduced  a 
data  model  that  states  the  dependency  of  the  clutter  reflec¬ 
tions  on  the  transmitted  signal.  Then,  we  developed  a  new 
polarimetric  detector  that  is  robust  against  space  and  time 
variations  of  the  clutter.  The  proposed  detector  decides  the 
presence  of  a  target  based  only  on  the  data  recorded  from 
the  cell  under  test,  without  requiring  secondary  data  or  prior 
knowledge  about  the  clutter.  We  studied  its  performance 
and  demonstrated  that  it  supports  the  CFAR  property.  We 
tested  this  detector  with  synthetic  and  real  data,  showing 
significant  performance  improvement  with  respect  to  other 
detectors.  An  additional  advantage  of  our  detection  algo¬ 
rithm  is  that  it  is  computationally  less  intense  because  it 
does  not  require  recording  and  processing  of  training  data. 


5.  APPENDIX 


The  target-to-clutter  ratio  and  clutter-to-noise  ratio  are 


TCR  = 


Ml 

tr(E)’ 


CNR  = 


tr(E) 


where  II  •  1 1  is  the  norm  of  the  vector. 
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